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Abstract. Description of parton interaction with the hot and dense nuclear matter created in 
f-O . heavy-ion collisions at high energies is a complex task, which requires a detailed knowledge of jet 

production in p+p and d+Au collisions. Measurements in these collision systems are therefore 
essential to disentangle initial state nuclear effects from cold nuclear matter effects, medium- 
induced kr broadening and jet quenching. To obtain complete description of the p+p (d+Au) 
collision it is also important to study particle production in the underlying event. The measured 
properties of underlying event can be used to tune the QCD based Monte-Carlo models. 
i— { • In this paper some of the recent results on jet and underlying event properties in p+p and 

d+Au collisions at ^snn = 200 GeV measured by the STAR experiment are presented. In 
particular, the preliminary results on charged and strange particle fragmentation functions in 
p+p collisions are discussed and confronted with PYTHIA simulations. Next, the measurement 
of inclusive jet spectrum and di-jet correlations in d+Au collisions is presented and compared 
with the measurements in p+p collisions to estimate the size of cold nuclear matter effects. 
Finally, the study of underlying event properties in p+p collisions is shown and compared to 
qq \ PYTHIA simulation. 
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1. Introduction 

Full jet reconstruction is considered as a promising tool for the quantitative study of properties 
of the hot and dense medium produced in heavy-ion collisions at RHIC energy. The complexity 
of heavy- ion collisions, where modifications of the fragmentation functions due to interaction 
of partons with the hot and dense medium are expected, requires inevitably a detailed 
understanding of jet production in p+p and d+Au collisions. Therefore measurements in these 
collision systems can be used as a baseline to disentangle initial state nuclear effects from cold 
nuclear matter effects, medium-induced fey broadening and jet quenching. In addition, it is 
also important to understand how the beam-beam remnants, multiple parton interactions and 
initial- and final-state radiation combine to produce particles in the underlying event (UE). 
The extracted underlying event properties serve as in important input for tuning of the QCD 
Monte-Carlo models such as PYTHIA. 

The first results on inclusive jet spectra in p+p collisions at y^ = 200 GeV measured by the 
STAR experiment at RHIC in 2003-2004 were published in [lj. The measured inclusive jet cross- 
section is in a very good agreement with NLO pQCD calculations over seven orders of magnitude. 
During the data taking in 2006 and 2008, the STAR experiment collected large data samples 
of p+p and d+Au collisions at y/sJm = 200 GeV. In this paper, the study of p+p collisions is 
extended by investigating fragmentation functions of charged and strange particles (Kg, A, A) 



as well as underlying event properties. Further, the paper discusses recent preliminary results 
on inclusive jet spectra and di-jet correlations in d+Au collisions. 

2. Dataset and analysis 

The preliminary results presented in this paper are based on analysis of p+p collisions at 
yfs =200 GeV measured in 2006 and p+p and d+Au collisions at ^Jsmn =200 GeV recorded 
in 2008 by the STAR experiment. To increase yield of hard processes, two special triggers 
utilizing information from Barrel Electromagnetic Calorimeter (BEMC) were used in addition 
to minimum bias trigger. In 2006, a "jet-patch" trigger defined as a A?y x A(f> = 1x1 patch 
of BEMC containing transverse energy Et > 8 GeV was used. In 2008, an online BEMC high 
tower (HT) trigger was used. To satisfy the HT trigger condition, at least one tower of BEMC 
must contain Et above a predefined threshold, which in case of the data presented here was 
4.3 GeV. In d+Au collisions, the Beam Beam Counter detector, located in the fragmentation 
region of Au nucleus, was used to select the 20% highest multiplicity events. 

Jets were reconstructed using sequential recombination algorithms kt and anti-kt as well 
as seedless infrared safe cone algorithm SISCone, which are implemented in the FastJet 
package [21 El H]. For the measurement of the neutral component of jets BEMC was used. 
The charged component of jets was measured by the Time Projection Chamber (TPC). Both 
TPC and BEMC have a full azimuthal coverage and their combined pseudorapidity acceptance 
is \r]\ < 0.95. This limits the fiducial acceptance of jets in pseudorapidity to 1-R, where R is the 
resolution parameter used in jet reconstruction. To avoid double counting of energy (electrons, 
minimum ionizing particles and possible hadronic showers in BEMC), the so called 100% 
hadronic correction was applied: the momentum of a charged particle from TPC was subtracted 
off the neutral tower energy in BEMC, if the measured charged particle track pointed to a BEMC 
tower. In order to minimize biases due to dead areas and possible background in BEMC, the 
neutral energy fraction (NEF) of the total jet energy was required to be within 0.1 < NEF < 0.9. 
For charged particle tracks, an upper transverse momentum cut-off pt < 15 GeV/c was applied 
due to uncertainty in TPC tracking performance at large transverse momentum. Detailed 
investigation of this limitation is in progress. 

In parallel to experimental data analysis, the analysis of jets and underlying event properties 
was also performed on simulated data. The simulations are based on PYTHIA 6.410, tuned to 
the CDF data at yfs =1.96 TeV (Tune A) and combined with a full STAR detector simulation 
based on GEANT. The simulated data were then analyzed using the same analysis chain as that 
for the measured data. We distinguish several levels of the simulated data: jet reconstruction at 
Monte Carlo particle level (PyMC) and at detector level including GEANT simulation (PyGe). 
In addition, to determine influence of background and its fluctuations on jet production in d+Au 
collisions, PYTHIA events were embedded in 0-20% highest multiplicity minimum bias d+Au 
events. This procedure was performed at detector level using reconstructed TPC tracks and 
signal in BEMC towers. This level of simulation is labeled as "PyBg". 

3. Fragmentation functions in p+p collisions 

Figure [Tj shows the uncorrected fragmentation functions of charged particles for jets with 
transverse momentum 20 < pt < 30 GeV/c measured in p+p collisions at yfs = 200 GeV using 
the jet-patch triggered data [S1E1E]- As this trigger creates a neutral energy fragmentation bias 
for the triggered jet, the fragmentation functions are shown only for the di-jet partner on away 
side. The fragmentation functions are plotted as a function of £ = ln(p^ /p^ adron ). The jets 
were reconstructed by three different jet algorithms (kt, anti-&£ and SISCone) with a resolution 
parameter R = 0.4 and R = 0.7, respectively. All charged tracks measured by TPC and neutral 
towers from BEMC with pt(Et) > 0.2 GeV/c were included in jet reconstruction. No significant 
dependence on jet algorithm is observed. The measured fragmentation functions are compared 
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Figure 1. Charged particle fragmentation function as a function of £ at detector level for 
jets with 20< px < 30 GeV/c measured in p+p collisions at yfs = 200 GeV. The jets were 
reconstructed by three different jet algorithms (triangles: k t , circles: anti-fct, squares: SISCone) 
with a resolution parameter R = 0.4 (left) and R = 0.7 (right) 0E1E]- The data are compared 
to PYTHIA simulations (histogram). 



to PYTHIA simulations described above. For both resolution parameters studied there is a good 
agreement between the data and simulation. This suggests that the NLO contributions which 
are not included in PYTHIA are small at RHIC. 

Further it is important to study in detail particle dependence of fragmentation functions 
as very little is known about identified fragmentation functions both experimentally and 
theoretically. While PYTHIA and NLO calculations agree reasonably well with the data on pion 
and proton production [8], theoretical description of strange particle production (K, A, 3) is a 
challenging task. The leading order calculations significantly underpredict the measured strange 
particle yields at intermediate and large transverse momenta [9]. Higher-order corrections 
to leading-order pQCD cross sections, which are in general large, can be accounted for by 
introducing a if- fact or, which renormalizes the leading-order pQCD cross section. It was found 
that a reasonable agreement of strange baryon particle production with PYTHIA at RHIC 
energy can be achieved by introducing K = 3 [9]. However, such a large value of K breaks the 
agreement of pion and proton data with PYTHIA. 

The STAR experiment has excellent capability to reconstruct strange particles such as Kq 
and A via their typical V0 decay channels with a very good signal-to- noise ratio at px > 1 GeV/c. 
Figure [2] shows the measured fragmentation functions of K s , A and A particles in p+p collisions 
at y/s = 200 GeV at detector level |10j . The fragmentation functions were calculated for three 
different jet px ranges: 10-15, 15-20, 20-40 GeV/c. The remaining small residual background 
under the mass peaks is not yet corrected for. The data are compared to polynomial fits to 
PYTHIA predictions. PYTHIA describes well the fragmentation functions of K s , which is in 
agreement with measurements of K fragmentation functions in e + + e~ collisions for similar jet 
energy range. For A and A fragmenation functions, PYTHIA predicts well the over-all yields, 
but it overestimates the yields at low £ and underestimates the yields at intermediate values of £. 
The integrated A/A and A-/K® ratios fro pr > 1 GeV/c in the studied jet px range are consistent 
with the values obtained from the minimum bias inclusive spectra when the same pt selection 
is applied [10J. This observation indicates that hard processes significantly contribute to K%, A, 
and A production at pt > 1 GeV/c. Further studies are needed to verify this. Corrections of 
the data to particle level as well as removal of remaining combinatorial background under mass 




1 5 < Reco Jet p T < 20 GeV/c 



20 < Reco Jet p T < 40 GeV/c 



[ STAR Preliminary 


" /+ +\ 


- * 

■ , jf*. i .... i ... 1^ ... . 



' STAR Preliminary 


^ A 


+- 


- ! ~*r 


: ,s*Ht--±-'. 

■ , ,4m , , , , i , , , , ri^_ 



5 log(jp T /fp T ) 



5 iog(Jp/fp T ) 



3 4 

5 log(jP/P T ) 



Figure 2. Fragmentation functions of Kg (squares), A (triangles) and A (upside-down triangles) 
as a function of the reconstructed jet px in p+p collisions at yfs = 200 GeV [ID]. The curves 
show polynomial fits to PYTHIA 6.4 predictions. Both data and simulations are shown at 
detector level. 



peaks are ongoing. 



4. Inclusive jet spectrum in d+Au collisions 

The inclusive transverse momentum spectrum of jets was measured in 10M 0-20% highest 
multiplicity d+Au collisions at ^snn = 200 GeV. All charged tracks measured by TPC and 
neutral towers from BEMC with px > 0.2 GeV/c were included in jet reconstruction. The jets 
were reconstructed with the anti-fcj jet finding algorithm with a resolution parameter R = 0.4. 
For subtraction of d+Au underlying event background, a method based on active jet areas [TT|[1"2"] 
was applied on event-by-event basis. The inherent asymmetry of the d+Au collision system 
results in r) asymmetric underlying background. The r\ dependence of underlying background 
was parametrized by a linear function and included in the background subtraction procedure. 

To correct for detector effects in d+Au collisions, the PYTHIA simulation of p+p collisions 
was used. In addition, a correction accounting for the lower TPC reconstruction efficiency in 
d+Au relative to p+p collisions was applied to avoid bias due to different Jet Energy Scale 
(JES). The TPC reconstruction efficiency in d+Au was determined from a single charged pion 
embedding into measured minimum bias d+Au events on detector level. The higher tracking 
efficiency in the PYTHIA simulation was then artificially lowered to match the value in d+Au 
collisions. This procedure was applied prior to jet finding at PyGe and PyBg level of simulations. 

To correct the jet spectrum to particle level, a bin-by-bin correction method was used. This 
correction procedure is based on the generalized efficiency determined as the ratio of PyMC 
to PyBg jet pt spectra, which is applied to the measured jet pt spectrum in d+Au collisions. 
The validity of this correction procedure relies on the fact that the shapes of the PyBg and the 
measured jet pt spectra are consistent. In [13J it was demonstrated that this condition is within 
statistical errors satisfied. 

The fully corrected inclusive jet px spectrum in d+Au collisions at y/SNN = 200 GeV is 
shown in Figure [3] [13, 14J. The systematic uncertainties are dominated by the JES uncertainty 
consisting of 10% systematic uncertainty in TPC charged particle reconstruction efficiency and 
5% in BEMC calibration uncertainty, added in quadrature. Ongoing studies with a full jet 
embedding into measured d+Au events at raw detector level will enable to decrease the JES 
uncertainty. 
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Figure 3. Inclusive jet px spectrum in 0-20% highest multiplicity d+Au collisions at 
y/SNN =200 GeV compared to {N^ in ) scaled p+p spectrum from [I]. The systematic errors 
(grey boxes and a dashed band) are explained in the legend. 



The jet spectrum in d+Au collisions is compared with the jet pt spectrum in p+p collisions 
from 2003-2004 for which a Mid Point Cone (MPC) jet algorithm with a cone radius R = 0.4 was 
used m . The effect of a slightly different r] acceptance on jet pt spectrum in case of p+p data 
was found to be less than 10% in the covered pt range. For comparison with the jet spectrum 
in d+Au collisions, the jet spectrum in p+p collisions was scaled by an average number of 
binary collisions {Nf,i n ) = 14.6+1.7 calculated from the Monte-Carlo Glauber model for 0-20% 
highest multiplicity d+Au events. The jet spectra in both collision systems are consistent within 
current systematic errors dominated by the JES and MC Glauber model uncertainties. A future 
comparison of the inclusive jet spectrum in d+Au collisions with that in p+p collisions measured 
during the same data taking period in 2008 is expected to further decrease the systematic 
uncertainties because they will be largely correlated. 



5. Nuclear k? broadening 

Azimuthal correlations of di-jets can be used to quantify size of nuclear kx broadening in d+Au 
collisions relative to p+p collisions. The total transverse momentum of a di-jet (kT,raw) can be 
determined from the di-jet opening angle A(f>: kj- = px sin(A(/>). To increase the di-jet yield, the 
HT triggered dataset with E+ > 4.3 GeV was used for both p+p and d+Au data measured in 
2008. Jets were reconstructed with the anti-kt algorithm using the resolution parameter R = 0.5. 
In addition, a lower transverse momentum cut pt > 0.5 GeV/c was applied to charged tracks 
and neutral towers to reduce background. For di-jet analysis the two most energetic jets in a 
given event (pt,i > Pt,2) satisfying 10 < pr,2 < 20 GeV/c were selected. To obtain a clean di-jet 
sample a cut on di-jet separation in azimuth \A<fi — ir\ < tt/3 was applied and measured kT,raw 
distributions were constructed as kr,raw = Pt ,i sin(A(/>) . 

To quantify the size of detector and background effects on kx,raw distributions, a di-jet 
analysis was performed on PYTHIA simulated data which is demonstrated in Figure H] together 
with a Gaussian fit of the simulated distributions |13| . 114 ^ fT5] . Comparison of PyMC, PyGe and 
PyBg distributions shows that the extracted Gaussian widths of kT,raw distributions are within 
statistical errors consistent. This suggests that detector and d+Au background effects in data 
can in first order be neglected. It should be however noted, that in case of PYTHIA simulations 
without background (PyMC and PyGe), the Gaussian fit does not describe the data well. 
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Figure 4. Distributions of kr,raw = Pt,i sin(A</>) for 10 < pr,2 < 20 GeV/c in MC simulated 
data PyMC (left), PyGe (middle) and PyBG (right). Jets were reconstructed with the anti-kt 
algorithm with a resolution parameter R = 0.5. The distributions were fit with a Gaussian fit 
function. 
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Figure 5. Measured distributions of fcr,™™ = Pt,i sin(Ac/>) for 10 < pt,2 < 20 GeV/c in p+p 
(left) and 0-20% highest multiplicity d+Au collisions (right) at ^/snn = 200 GeV. Jets were 
reconstructed with the anti-kt algorithm with a resolution parameter R = 0.5. The distributions 
were fit with a Gaussian fit function. 
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The fcr.raw distributions for data from p+p and d+Au collisions are shown in Figure The 

2.8+0.1 (stat.) GeV/c in p+p 

3.0+0.1 (stat.) GeV/c in d+Au collisions p31 [M [[5] . The widths in 
both collision systems are consistent within statistical errors indicating very small nuclear kx 
broadening. In case of p+p collisions, the Gaussian fit does not describe the shape of the A>r,raw 
distribution well similarly as observed for the PyMC and PyGe simulated data. The precise 
shape of the A>r,raw distribution is currently under investigation. Varying the size of |A</> — vr| cut 
for back-to-back di-jet selection between 0.5 and 1.0, the systematic uncertainty on the extracted 
Gaussian widths of A>r,raw distribution was estimated to be 0.2 GeV/c. The effects of pt,2 range 
and jet finding algorithms were found to be negligible in both studied collision systems. 




Figure 6. Cartoon of the event areas ("toward", "away", "transverse") in azimuthal angle 
difference relative to the leading jet (A</>). 



6. Underlying event studies in p+p collisions 

A full description of particle production in p+p collisions is a complex task which requires to 
incorporate both production of jets as well as production of particles in the so called underlying 
event within the same model. The underlying event (UE) is in general defined as everything but 
the hard scattering. It has contributions from soft and semi-hard multiple parton interactions 
(MPI), initial and final state radiation effects (ISR, FSR) and beam-beam remnants. As the 
study of UE presented below is focused on mid-rapidity, the contribution of beam-beam remnants 
to the UE properties can be neglected. 

The UE study presented here follows the analysis technique developed by the CDF 
Collaboration at Tevatron |16} 117] . After reconstructing jets in a given event, the event 
is split into four sections schematically depicted in Figure [6) The sections are defined by 
their azimuthal angle difference with respect to the leading jet axis (A0). The "toward" 
region is the region within \A<p\ < 60° around the leading jet axis and the "away" region 
is located at \A(f>\ > 120°. The remaining two regions labeled as "transverse" are placed at 
60° < | A(f>\ < 120° relative to the leading jet axis. The underlying event contributes to all of the 
mentioned regions, but to access its properties the two "transverse" regions are used as there the 
contamination from jets is assumed to be small. According to the charged particle multiplicity 
(N c h) the two transverse regions are further labeled as "TransMax" and "TransMin", where 
A c /i (TransMax) > N c h (TransMin) in a given event. This selection enhances contribution of 
initial and final state radiation components of the hard scattered parton in TransMax region 
relative to TransMin region. The sensitivity to various physics processess in UE can be achieved 
by comparing properties of UE for a leading jet sample and a sample of events containing di- 
jets, which suppresses the probability of ISR and FSR effects. For this study events with two 
reconstructed jets with p^" ay /pi^ ad > 0.7 and |A(/>| > 150° were selected. 

Figure [7] shows the uncorrected charged particle density and mean transverse momentum in 
TransMin, TransMax and away side jet regions as a function of leading jet pt in the di-jet sample. 
Jets were reconstructed with the kt algorithm with the resolution parameter R = 0.7 [51E1E]- 
No significant dependence on jet algorithm was observed. The charged particle density as well 
as {pt) are approximately independent of leading jet pt, while the away-side jet charged particle 
density and [pt) are rising with jet pt, as expected. Compared to data measured at Tevatron 
(^/s = 1.96 TeV), the values of chaged particle density and (pr) in the underlying event at RHIC 
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Figure 7. The uncorrected charged particle density (left) and charged particle mean p? (right) 
in the "di-jet" data set for the away side jet, TransMin and TransMax regions as a function of 
reconstructed leading jet pr- The jets were reconstructed with the kt algorithm with a resolution 
parameter R = 0.7. The dashed histogram indicates the prediction from PYTHIA simulations 
(Tune A, e = 0.25). The data and simulation are presented at detector level. 



energy are in general lower. 

The measured data in Figure [7] are compared to PYTHIA simulation (Tune A) with MPI 
scaling factor e = 0.25 based on recent measurements at Tevatron energies {yfs = 630, 1800 
and 1960 GeV) |18| . PYTHIA describes the main features of the data well, however a closer 
investigation reveals that PYTHIA slightly underestimates the charged particle multiplicity and 
(pt) in TransMax and TransMin regions while overestimating these quantities within the jet. 

In Figure [8] the measured charged particle density in UE is investigated closer by comparing 
the results for the transverse regions between the leading and di-jet samples [HI El [7]. No 
significant difference between the jet samples is observed at RHIC energy which points to the fact 
that the hard scattered partons do not emit significant amount of large angle ISR and FSR. This 
observation is very different to that at Tevatron, where in p + p collisions at yfs = 1.96 TeV the 
ratio of charged particle density in leading/di-jet samples was found to be about 0.65 [TTJ. The 
data in Figure [8] are compared with an expected charged particle density in case the multiplicity 
in the UE would follow a simple Poisson distribution with an average value of 0.36. This simple 
simulation suggests that at RHIC energy, the splitting of the TransMax and TransMin values is 
mostly due to the sampling of the total charged particle multiplicity measured. 



7. Summary 

In summary, various properties of jets in p+p and d+Au collisions at y/SNN = 200 GeV were 
discussed and compared to PYTHIA simulations. The measured fragmentation functions of 
charged particles and identified strange particles at detector level show a good agreement with 
PYTHIA simulations for charged and Kg particles, while for A and A PYTHIA fails to describe 
details of the fragmentation pattern. 

The measured inclusive jet pt spectra and nuclear kx broadening in d+Au collisions from 
2008 do not show any significant cold nuclear matter effects when compared to the measurements 
in p + p collisions within current statistical and systematic uncertainties. These uncertainties 
will be further decreased by comparing the data to a new jet cross section measurement in 
p+p collisions from 2008 and improvements in the tracking efficiency determination from jet 
embedding in real d+Au data. 
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Figure 8. The uncorrected charged particle density in the TransMin and TransMax regions as 
a function of reconstructed leading jet pt using SIS Cone algorithm with a resolution parameter 
R = 0.7 shown separately for a leading jet and back-to-back jet samples. The data and simulation 
are presented at detector level. The solid line indicate expected charged particle density in case 
of simple Poisson distribution of charged particles in UE. 



First studies of underlying event properties in p+p collisions show that its charged particle 
density and mean p? are largely independent of jet transverse momentum and the contributions 
from initial and final state radiation are small. Studies of UE in d+Au collisions are underway 
and are expected to bring further insights into details of particle production mechanism in cold 
nuclear matter. 



Acknowledgements 

This work was supported in part by grants LC07048 and LA09013 of the Ministry of Education 
of the Czech Republic. 



References 

[1] B. Abelev et al. (STAR Collaboration), Phys. Rev. Lett. 97 (2006) 252001. 

[2] M. Cacciari, G. Salam, Phys. Lett. B 641 (2006) 57. 

[3] G. P. Salam, G. Soyez, JHEP 0705 (2007) 086. 

[4] M. Cacciari, G. Salam, G. Soyez, JHEP 0804 (2008) 063. 

[5] H. Caines et al. (STAR Collaboration), Nucl. Phys. A 830 (2009) 263C. 

[6] H. Caines et al. (STAR Collaboration), arXiv:0910.5203l 

[7] H. Caines et al. (STAR Collaboration), Nucl. Phys. A 855 (2011) 376. 

[8] J. Adams et al. (STAR Collaboration), Phys. Lett. B 637 (2006) 161. 

[9] B. I. Abelev et al. (STAR Collaboration), Phys. Rev. C 75 (2007) 064901. 

[10] A. Timmins et al. (STAR Collaboration), J. Phys. G 37 (2010) 094035. 

[11] M. Cacciari, G. Salam, Phys. Lett. B 659 (2008) 119. 

[12] M. Cacciari, G. Salam, G. Soyez, JHEP 0804 (2008) 005. 

[13] J. Kapitan et al. (STAR Collaboration), Nucl. Phys. A 855 (2011) 412. 

[14] J. Kapitan et al. (STAR Collaboration), J. Phys. Conf. Ser. 270 (2011) 012015. 

[15] J. Kapitan et al. (STAR Collaborati on), Nucl. Phys. A 830 (2009) 619C. 

[16] R. Field et al. (CDF Collaboration), |hep-ph/0510198| 

[17] A. Cruz (CDF Collaboration), PhD Thesis, University of Florida, USA (2005). 

[18] R. Field, AIP Conf. Proc. 828 (2006) 163. 



